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Abstract 

Recent improvements to the Cascade-Exciton Model (CEM) of nuclear reactions are briefly de- 
scribed. They concern mainly the cascade stage of reactions and a better description of nuclei during 
—i | the preequilibrium and evaporation stages of reactions. The development of the CEM concerning fis- 

q ' sion is given in a separate talk at this conference. The increased accuracy and predictive power of the 

CEM are shown by several examples. Possible further improvements to the CEM and other models are 
discussed. 

^- : 

1. Introduction 



During the last two decades, several versions of the Cascade-Exciton Model (CEM) [JTp of nuclear 
reactions have been developed at JINR, Dubna. A large variety of experimental data on reactions induced 
by nucleons 0], pions 0], and photons [Q] has been analyzed in the framework of the CEM and the general 
validity of this approach has been confirmed. Recently, the CEM has been extended by taking into account 
the competition between particle emission and fission at the compound nucleus stage and a more realistic 
calculation of nuclear level density []50. In the last few years, the CEM code has been modified [j^, [70 to 



X 



calculate hadron-induced spallation and used to study Qq]-QT20 about 1000 reactions induced by nucleons 
from 10 MeV to 5 GeV on nuclei from Carbon to Uranium. 

At present, several versions of the CEM code are used at a number of national laboratories and uni- 
versities to solve different problems. Several versions of the CEM code have been implemented with or 
without modifications into a number of well-known transport codes, e.g., in: MARS fll3tl, HETC [ [14| , [15]], 
TIERCE [|Icj{], MCNPX [[T70, and is currently being incorporated into the LLNL radiation-transport code 
COG flITI[]. The preequilibrium part of the CEM has been implemented either as the initial Modified Exci- 
ton Model (MEM) code MODEX [ |T9[ ] or as the later version from 02Tj| ] wholly or with some modifications 
into the tranport codes HETC96 [[21} [Bj], SHIELD 0, HETC-3STEP @ [TJ], HETC-FRG [TJ], 



CASCADE [250, and into the intranuclear cascade-preequilibrium-evaporation-fission code INUCL [260. 



The recent International Code Comparison for Intermediate Energy Nuclear Data [ [270 has shown that 
the CEM adequately describes nuclear reactions at intermediate energies and has one of the best predictive 
powers for double differential cross sections of secondary nucleons as compared to other available models 
(see Tabs. 5 and 6 in the Report [^70 and Fig. 7 in Ref. fl28[]). As an example, in Fig. 1 we show experi- 
mental spectra of neutrons emitted from interactions of 1.5 GeV protons with C, Al, Fe, In, and Pb [[290 



compared with calculations using a version of the CEM as realized in the code CEM95 flZOQ. We see that 
CEM95 reproduces the data very well. In effect, these results can be regarded as a prediction of CEM95, 
as the experimental data were published in 1997 while the calculations are made without any changes 
of the code and using the previously determined "preferred" set of input parameters (RM=1.5, IFAM=9, 
IDEL=0, ISH=6, ISHA=1, IDELTA=1; see notations and details in the CEM95 code manual 02^]); these 
were fixed in 1995. 
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Fig. 1. Comparison of measured [^] double differential cross sections of neutrons from 1.5 GeV protons on C, 
Al, Fe, In, and Pb with CEM95 calculations. 



Unfortunately, not all CEM95 results agree so remarkably well with measuremetnts. As an example, 
Fig. 2 shows also experimental spectra of neutrons emitted from interactions of 256 MeV protons with 
Al [Q, Fe and Zr Q, and 1.5 GeV tt+ with Fe [|32|]. Calculations with CEM95 using the same 
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Fig. 2. Comparison of measured neutron spectra from 256 MeV protons on Al Q30[] , Fe J30|], and Zr [J31[], and 1.5 
GeV 7T + on Fe [^] with CEM95 calculations (solid histograms) and with results of LAHET [33] for Al from [30], 
prediction of the recently improved version of the Cugnon et al. INC for Fe, HETC-KFAII 0] calculations 



for Zr from Q31[], and a Japanese version [ |36| ] of HETC [ ]15[ ] using modified NN cross sections to take account of the 
in-medium effects for tt + (1.5 GeV)+Fe from Ref. [32], shown by dashed histograms. 



fixed set of input parameters as for Fig. I are shown by solid histograms. For comparison, results obtained 



with other well-known codes are shown as well: LAHET [ |33| ] calculations for Al from Q30|], a recent 
improved version of the Liege INC model by Cugnon et al. [ J3"3[ ] for Fe, a Jiilich version of HETC [ |T"5| ] 
as realized in the code HETC-KFAII j35t | for Zr from [^Tj], and a recent Japanese version [[36]] of HETC 
using modified NN cross sections to take account of the in-medium effects for 7r + (1.5 GeV)+Fe^n+... 



from Ref. [f32p. 

One can see that for these reactions the agreement of CEM95 results with the measured neutron spectra 
is not as good as that shown in Fig. 1; in places the difference is bigger than a factor of 2. This seems to 
be a characteristic of all the models compared here. 

For any model, to predict yields of isotopes produced at intermediate energies is much more difficult 
than to calculate spectra of emitted particles []37|]. As one can see from Figs. 3 and 4 (adapted from [[381]), 
CEM95 describes the majority of yields for isotopes produced in the spallation region quite well and no 
worse than other well-known codes, although for some individual isotopes the discrepancy with the 
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Products in Bi-209 irradiated with 0.13GeV protons 




Fig. 3. Comparison between experimental data 
LAHET 
protons. 



(filled circles) and calculations with the codes CEM95 [2C], 
3j, and INUCL [M (open circles) of yields of the indicated isotopes from 209 Bi irradiated with 130 MeV 
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Fig. 4. The same as Fig. 3 using instead calculations with the CASCADE [25], HETC [15], HMS-ALICE [39] 



codes, and YIELDX phenomenological systematics [40]. 
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cross sections, measurements is quite large. 

The recent projects for Accelerator Transmutation of Waste (ATW) and Accelerator Production of 
Tritium (APT) have revived interest in accurate nuclear reaction data at intermediate energies (see, e.g., 
pl]]). Not all needed data can be measured, therefore reliable models are required to provide the necessary 
This is our motivation for attempting to improve the CEM into a model capable of predicting reliable 
nuclear cross sections for arbitrary targets in a wide range of incident energies. 

The philosophy of our work is not just to include a few more parameters, fitting them to describe the 
available data, a commonly used practice in development of codes for applications. Instead, we are fur- 
ther developing the CEM by progressively incorporating features of previously neglected physics, striving 
meanwhile not to destroy the current broad strengths of CEM95 and its good predictive power. To in- 
vestigate the dependence of theoretical results on the physics incorporated in CEM95 and to identify the 
improvements to the CEM and to other similar models which are of highest priority, we have performed a 
detailed analysis of more than 600 excitation functions for proton-induced reactions on 19 targets from C 
to Au at energies from 10 MeV to 5 GeV []7]]. After this analysis, we created a list of potential improve- 
ments to the CEM95 code [Q]. Some these have already been incorparated and we describe them briefly 
in the present paper. Nevertheless, as we have not yet completed our work, some of the results presented 
below are still preliminarily. 

2. The Main Concepts of the Model 

A detailed description of the CEM may be found in Ref. [JTp and of its extended version, as realized in 
the code CEM95, in Refs. [§, 0]; therefore, we mention here only its basic assumptions and features modi- 
fied in the present work, for the sake of clarity. The CEM assumes that reactions occur in three stages. The 
first stage is the intranuclear cascade in which primary and secondary particles can be rescattered several 
times prior to absorption by, or escape from the nucleus. The cascade stage of the interaction is described 
by the standard version of the Dubna intranuclear cascade model (INC) [|42|]. The excited residual nucleus 
remaining after the emission of the cascade particles determines the particle-hole configuration that is the 
starting point for the second, preequilibrium stage of the reaction. The subsequent relaxation of the nuclear 
excitation is treated by an extension of the Modified Exciton model (MEM) [[T^] of preequilibrium decay 
which also includes a description of the equilibrium evaporative third stage of the reaction. 

All the cascade calculations are carried out in a three-dimensional geometry. The nuclear matter den- 
sity p(r) is described by a Fermi distribution with two parameters taken from the analysis of electron- 
nucleus scattering, namely 

p(r) = p p (r) + p n (r) = p {l + exp[(r - c)/a]} , (1) 

where c = 1.07A 1 / 3 fm, A is the mass number of the target, and a = 0.545 fm. For simplicity, the target 
nucleus is divided by concentric spheres into seven zones in which the nuclear density is considered to be 
constant. The energy spectrum of the target nucleons is estimated in the perfect Fermi gas approximation 
with the local Fermi energy T F (r) = h 2 [3n 2 p(r)] 2 / 3 /(2m N ), where ui]\f is the nucleon mass. The influence 
of intranuclear nucleons on the incoming projectile is taken into account by adding to its laboratory kinetic 
energy an effective real potential V, as well as by considering the Pauli principle which forbids a number 
of intranuclear collisions and effectively increases the mean free path of cascade particles inside the target. 
For incident nucleons V = Vjv(r) = T F {r) + e, where T F {r) is the corresponding Fermi energy and e 
is the mean binding energy of the nucleons (e ~ 7 MeV [|42|]). For pions, in the Dubna INC one usually 
uses a square-well nuclear potential with the depth V w — 25 MeV, independently of the nucleus 
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and pion energy. The interaction of the incident particle with the nucleus is approximated as a series of 
successive quasifree collisions of the fast cascade particles (N or n) with intranuclear nucleons: 

NN -> NN, NN -> irNN, NN -> TTi, • • ■ , t^AW , 

ttN^ttN, irN -> ir 1 ,---,7r i N (i > 2) . (2) 

To describe these elementary collisions, we use experimental cross sections for the free iViV and 7riV 
interactions, simulating angular and momentum distributions of secondary particles by special polynomial 
expressions with energy-dependent coefficients, and take into account the Pauli principle. 

Besides the elementary processes (2), the Dubna INC also takes into account pion absorption on nu- 
cleon pairs 

7T + [NN] -> NN. (3) 

The momenta of two nucleons participating in the absorption are chosen randomly from the Fermi distri- 
bution, and the pion energy is distributed equally between these nucleons in the center-of-mass system of 
the pion and nucleons participating in the absorption. The direction of motion of the resultant nucleons in 
this system is taken as isotropically distributed in space. The effective cross section for absorption (let us 
speak below, for concreteness, e.g., about 7r _ ) is estimated from the experimental cross-section for pion 
absorption by deuterons 

a(n~ + u np" —> nn) = W ■ a(n~ + d — > nn) . (4) 

The quantity W can depend on the pion energy T w , on the characteristics of the target nucleus, the point 
where the pion is absorbed, and on the spin-isospin states of absorbing pairs (see details and references 
in @D. 

3. Study Results 

A. Elementary Cross-Sections. In the Dubna INC model [|42|] used in CEM95, the cross sections for 
the free NN and tcN interactions (2) are approximated using a special algorithm of interpolation/extrapolation 
through a number of picked points, mapping as well as possible the experimental data. This was done very 
accurately by the group of Prof. Barashenkov using all experimental data available at that time, about 30 
years ago. Currently, the experimental data on cross section is much richer, therefore we decided to revise 
the approximations of all elementary cross sections used in CEM95. We started with collecting all pub- 
lished experimental data from all available sources. Then, we developed an improved, as compared with 
the standard Dubna INC [@2j], algorithm for approximation of cross sections and developed simple and fast 
approximations for elementary cross sections which fit very well presently available experimental data not 
only to 5 GeV, the upper recommended energy for the present version of the CEM, but up to 50-100 GeV 
and higher, depending on availability of data. So far, we have such approximations for 34 different types 
of elementary cross sections induced by nucleons, pions, and gammas. Cross sections for other types of 
interactions taken into account in the CEM are calculated from isospin considerations using the former as 
input. 

We consider this part of our CEM improvement as an independently useful development, as our ap- 
proximations are reliable, fast, and easy to incorporate into any transport, INC, BUU, or Glauber-type 
model codes. For example, our new approximations recently have been successfully incorporated by 
N. Mokhov into the recent improved version of the MARS code system at Fermilab [JT3"|]. 
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Fig. 5. Energy dependence of the np, n + p, and 7^ total cross sections and of the np —>■ ppir , ir + p — > ir + pTr°, 
and 7£> — > mr + i: ones. Experimental points are from our compilation [43]. Solid lines are results of our present 



approximations; dashed lines show the standard Dubna INC approximations Q42[ ] used in CEM95. 
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We are going to publish our compilation of experimental data and our approximations of these cross 
sections in a separate paper and to make them available to users through the Web [|43|], as well as to 
make available to users our code for these approximations. An example of 6 compiled experimental cross 
sections together with our new approximations and the old approximations from CEM95 is shown in 
Fig. 5. We see that our new approximations describe indeed very well all data. Although presently we 
have much more data than 30 years ago when Barashenkov's group produced their approximations used 
in CEM95, for a number of interaction modes like the total cross sections shown in the left panel of Fig. 5, 
the original approximations also agree very well with presently available data, in the energy region where 



the Dubna INC was developed to work. This is a partial explanation of why the old Dubna INC p2fl and 
the younger CEM95 [ PD[ ] work so well for the majority of characteristics of nuclear reactions. 

On the other hand, for some modes of elementary interactions like the ones shown in the right panel 
of Fig. 5, the old approximations differ significantly from the present data, demonstrating the need for our 
present improvements for a better description of all modes of nuclear reactions. 

B. Effects of Refraction and Reflection. In CEM95, the kinetic energy of cascade particles is in- 
creased or decreased as they move from one potential region (zone) to another, but their directions remain 
unchanged. That is, in our calculations, refraction or reflection of cascade particles at potential boundaries 
is neglected. 

To understand how our results depend on effects of refractions and reflections we performed calcu- 
lations for a number of reactions with a modified version of CEM95 taking into account refractions and 
reflections as described in the monograph p2"| ] and realized by A. S. Iljinov [fBJ]. An example of our re- 
sults is shown on Fig. 6. From these and other similar results obtained for other reactions we conclude that 
refractions and reflections affect only weakly the spectra of secondary particles and somewhere more the 
total inelastic/elastic cross sections. As one could expect in advance, at low incident energies, refractions 
and reflections more strongly affect the calculated spectra but on the whole, CEM95 results without taking 
into account refractions and reflections agree better with the data, for the majority of analyzed reactions. 
Therefore, we prefer for our new version of the CEM to ignore refractions and reflections, as in the stan- 
dard Dubna INC. We note that similar results were obtained earlier by Chen et al. [07]]. 



C. Pion-Nucleus Potential V w . It is natural that theoretical characteristics of a nuclear reaction involv- 
ing pions depend on the pion-nucleus potential used in calculations. This problem has a history of almost 
half a century; it is not solved completely even now, and its review is beyond the scope of the present 
paper. 

As mentioned above, CEM95 uses a constant attractive pion-nucleus potential of 25 MeV. This value 
was obtained 30 years ago by Barashenkov, Gudima, and Toneev [ |48| ] from an analysis of available data 
and was suggested as a basis for the Dubna INC. At the same time, analyzing pion-nucleus reactions 
at energies below 300 MeV and nuclear absorption of stopped antiprotons with the Dubna INC, Iljinov, 
Nazaruk, and Chigrinov suggested [ |49| ] that one can achieve a better agreement with the data using V w = 0. 
Making the situation more confusing, in a recent extension of the Dubna INC for photonuclear reactions 
at energies up to 10 GeV by Iljinov et al. [|50|], a better description of pion photoproduction was achieved 
for K = 35 MeV. 

Since previous INC studies of the values of V n are based on analysis of quite scanty and old data, 
we decided to return once more to this question using the recent measurements at LAMPF of non-charge 
exchange (NCX) [|5T|], single charge exchange (SCX) [ |33| , [52p, and double charge exchange (DCX) [ |53"| ] 



pion production on targets from C to Bi in the enegy region from 120 to 500 MeV. 
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Fig. 6. Measured double differential spectra and angular distribution of ir° from 7r + (165 MeV)+ 16 [45] and 
Barashenkov's systematics p6| ] for the experimental 7r + +Ca total inelastic cross section compared with CEM95 
calculations without (solid) and with refractions and reflections (dashed). 



We perform calculations for a large number of NCX, SCX, and DCX reactions using several values 
for V n . An example of our results is shown in Fig. 7. Our analysis shows that total inelastic cross sections 
and mean multiplicities of secondary pions, as well as NCX and SCX pion spectra at incident energies 
well above the A(1232) resonance region depend only slightly on the value of V n . But at lower incident 
energies corresponding to the A (1232) region the shapes of all NCX, SCX, and DCX spectra change 
significantly with changes in the value of V w . Unfortunately, none of the fixed values of the V w we try 
can provide a good simultaneous description of all data. Nevertheless, the best overall agreement with the 
different data we analyse is achieved using the original value of V n = 25 MeV. 

Other recent studies show that V n is a complex function of pion energy T n , nuclear target (A, Z), radius 
vector r, and on the isospin of the pion. Ideally, for any INC model it would be useful to have a simple 
approximation for V n (A, Z, T n , r) that could be used in simulations. We do not know any simple system- 
atics for V w , and it remains still an input parameter. For the moment, we choose V w = 25 MeV, with the 
understanding that this is not a final solution to the problem. 

D. Pion Transparency. Transparency means an inhibition of the interaction of a pion with subsequent 
nucleons in the target nucleus for a certain time after creation of the pion. Recently, attention was called 
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Fig. 7. Comparison of experimental NCX [51 ] and SCX [45] pion spectra for interactions of 500 MeV tt + with 
C and 165 MeV tt + with O with standard CEM95 results (V n = 25 MeV; solid histograms) and with calculations 
using Vtt = MeV (dashed histograms). 



to this possible new physical effect after a study of the previously mentioned NCX pion production mea- 
surements at LAMPF [|51]] with the INC model by Gibbs et al. (GINC) @. GINC failed to reproduce 
the data, strongly underpredicting the measured yields at pion energies near 200 MeV. When pions in- 
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volved in (tt, 2n) reactions were not allowed to interact with the nucleus for a time equal to 2 fm/c, a 
considerable improvement in the agreement with data was achieved |PTj], This fact was interpreted as 
a new phenomenon of nuclear pion transparency with several possible exotic sources like: a need for a 
hadronization time for the wave function of pions involved in pion production to settle into an eigenstate, a 
mechanism which allows the pions to be transported through the nucleus by weakly interacting o mesons, 
or renormalized particle masses at central nuclear density [|T]]. Recently, we analyzed [ |55| ] SCX n° spec- 
tra from interactions of negative pions with 12 C and 7 Li at 467 and 500 MeV using the same GINC [|54|] 
and we got results very similar to the ones on NCX spectra published by Zumbro et al. [|51|]; we saw also 
an indication of some pion transparency. 

Therefore we decided to implement this new piece of physics in the CEM, if appropriate. The Dubna 
INC used in CEM95 does not take into account the time of intranuclear elementary interactions, instead, it 
considers the coordinates of these interactions. So, instead of a hadronization time, we can require that the 
mean free path of secondary pions produced in (tt, 2w) reactions to be longer than a fixed value, e.g., 2 fm. 
Including a transparency distance of 2 fm for pions from (tt, 2tt) interactions into the CEM improves the 
description of some of the NCX, SCX, and DCX pion spectra at pion incident energies around the A (1232) 
resonance region. An even better agreement with some of the data [|5l]]-[^3]] is achieved by imposing a 
transparency of 2 fm for pions created in a (it, 2ir) interaction and, in addition, a transparency of 1.2 fm 
for all other pions involved in a reaction. 

As we mentioned in the beginning, the aim of our work is to incorporate progressively in CEM95 new 
physical ingredients without destroying its present wholeness and good predictive power. Therefore we 
consider it necessary to analyze not only pion-induced NCX, SCX, and DCX reactions, but also other types 
of intermediate energy nuclear reactions where pion production, absorption, and emission are important 
channels. We study a large variety of characteristics of different types of nuclear reactions from a unique 
point of view, without fitting any parameters except for several versions of pion transparencies. Two results 
of this study, one "good" and one "bad", are shown in Fig. 8. The general conclusion we drawn from our 
analysis is that imposing some transparency for pions in CEM95 improves description of parts of the NCX, 
SCX, and DCX pion spectra around the A(1232) resonance region (see the upper plot in Fig. 8) as well as 
of the 100-200 MeV region of backward nucleon spectra from nucleon- and pion-induced reactions. At 
the same time, incorporation of such pion transparency destroys the self-consistency of the CEM95 and 
results in a clear overestimation of pion production at forward angles (see the lower plot in Fig. 8) for 
different reactions. We interpret these results as an indication that we need a further improvement of the 
treatment of pion-nucleus interactions in CEM95 (e.g., inclusion of A and other resonances explicitly as 
cascade participants, a better approximation of experimental tiN and NN angular distributions and of pion 
absorption in the CEM) rather than as evidence of new physics related with some exotic pion transparency. 



E. Pion Absorption. The pion absorption cross section on nuclear pairs is treated in CEM95 according 
to Eq. (4), using experimental data for absorption on deuterons. Calculations by many authors reveal that 
an overall satisfactory description of different experimental data can be obtained over a large range of 
pion energies and target nuclei, provided that the "effective" approximation W ~ const = 4 is used. The 
probability of absorption of pions on different pp, np, and nn pairs is taken to be the same, whenever the 
absorption is allowed by the charge conservation law. This is quite a rough approximation and we have 
encountered problems using it to study pion-nucleus interactions in the A resonance region [J3p. 

If we assume that pions are absorbed on nucleon pairs only through the production with subsequent 
absorption of A isobars, tt[NN] — ► AN — ► NN, than from the isospin decomposition one gets the 
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Fig. 8. Comparison of experimental SCX [45] it spectra from 400 MeV 7r on different targets [52] and of tt + 
spectra from 562.5 MeV neutrons on Cu J56]] with results of CEM95 calculations with pion transparency, i.e., with 



an imposed > 2 fm for pions from (it, 2ir) interactions and > 1.2 fm for all pions after the first 7riV interaction 
(solid histograms) and with standard CEM95 results (dashed histograms). For comparison, the thick dashed lines 
show results of LAHET (Bertini INC) [0] calculations from 
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following absorption ratios Q57Q: 



a^+inn) — > np a w o(np) — > np <j n o (nn) — > nn 

n ) ! = 0.083, - \ , = 0.440, - \ { = 0.140 , 

v<k+ {np) —>■ pp 0-^+ {np) — > pp 0V+ {np) — > pp 

= 0.140, ^44^ = 0-083, 5^ = 1- ( 5 ) 
0V+ {np) —> pp a w + {np) — > pp 0^+ (np) — > pp 

To test this approximation, we started with the compilation of all presently available experimental 
cross sections for pion absorption on deuterons. Then, we estimate a n -{np) — > nn using Eq. (4) with 
W = 4. All other absorption cross sections, for other types of nucleon pairs and/or charges of pions 
are calculated according to Eq. (5). With such modifications, we study several characteristics of nuclear 
reactions determined mainly by pion absorption on nucleon pairs. 

Pion absorption makes a large contribution to spectra of secondary nucleons emitted at very backward 
angles in the energy region around 140 MeV (see, e.g., [|2|]). Our first results, for backward proton spectra 
from proton-induced reactions, are very encouraging (see the upper plot in Fig. 9); using the ratios (5) for 
pion absorption significantly improves the agreement of calculated spectra with the data. But when we 
try these ratios for neutron spectra we get, as in the case of pion transparency, an opposite result; neutron 
yields from the pion absorption mode is strongly overestimated when using the ratios (5). We note that a 
variation of the absolute normalization of the pion absorption cross section, i.e., of W would not solve this 
problem, since if we use a smaller value for W to adjust the neutron yields, then the proton yields from 
pion absorption will be too low. This negative result leads us to reject the ratios (5) for CEM and to return 
to equal probabilities of pion absorption on different nuclear pairs. 

In Fig. 9, both proton and neutron spectra at backward angles calculated with CEM95 (i.e., using equal 
probabilities for pion absorption on different pairs) lie above the data by a factor of 2 at energies around 
100 MeV. These nucleons come in our model mainly from pion absorption. This means that pion absorp- 
tion cross sections used for these targets were too high and W in Eq. (4) should be smaller. Similar results 
were obtained with CEM95 for other light and medium targets. To improve this situation, we re-examined 
the value of W. Fortunately, we presently have reliable spectra of both protons and neutrons measured 
at backward angles for the same targets by the same experimentalists at energies around 100-200 MeV 



Q62Q. In Fig. 10, we show part of our results: experimental spectra [ p2| ] of neutrons and protons emitted 
at 140° from the interaction of 1 GeV protons with different nuclei are compared with standard CEM95 
calculations and using twice lower pion absorption cross sections, i.e., W = 2 for all targets. One can 
see that for medium and light nuclei, the results obtained using W = 2 agree much better with data as 
compared with the standard CEM95 results (W = 4). From these and similar results obtained for other 
targets at other intermediate energies we can conclude that W seems to depend on the atomic number of 
the target; increasing from about 2 for light targets like C to about 4 for heavy targets like Pb. Fig. 11, 
where we show results obtained with W = 2 for C, W = 3 for Fe and Ni, and W = 4 for Pb, serves 
as a confirmation of this conclusion: Introduction of a pion absorption cross section increasing with A 
on nuclear pairs significantly improves the agreement with data both for neutron and proton spectra at 
backward angles around 100-200 MeV. Similar results were obtained for other reactions. We have not yet 
fixed the value of W in our new version of the code, as we wish to perform several more test calculations 
to determine W more accurately. But we already see that we will need to use an A-dependent function for 
W, which increases from about 2 for C to about 4 for Pb. 
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Fig. 9. Experimental proton spectra from 500 MeV p + Ni [ pQ ] and neutron spectra from 800 MeV p + C [ pip 
compared with standard CEM95 results and with modified calculations of pion absorpion cross sections according 



to Q57Q (noted in this figure as CEM98) as described in the text. CEM97 results are obtained with the same pion 
abosrption cross sections as the CEM95 ones, but with new elementary NN and nN cross sections and new nuclear 
masses. 
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Fig. 10. Experimental spectra of neutrons and protons emitted at 140 degrees from interactions of 1.0 GeV 
protons with C, Al, Cu, Sn, and Pb j62| ] compared with standard CEM97 results (dashed histograms) and with cal- 
culations using one-half thepion absortpion cross sections, i.e., W = 2, or: a a bs('K+ [NN]) = 2a a b s ('K + d]) (noted 
here as CEM98). 



F. Other Enhancements, CEM97a Code, and Further Work. We have made a number of other 
improvements to CEM95 concerning the description of the preequilibrium and evaporation stages of reac- 
tions. This includes incorporating a complete experimental mass table extended by the Moller-Nix mass 
calculations [|58l], along with corresponding ground-state microscopic corrections and pairing energies. We 
have also implemented a new level-density approximation which incorporates the calculated microscopic 



corrections and pairing energies and angular momentum dependent macroscopic fission barriers [p9|], 
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Fig. 11. Comparison of experimental neutron (proton, for 500 MeV protons on Ni) spectra from interactions of 
protons of 800 MeV on C @, 500 MeV on Ni [§(}, 1.5 GeV on Fe and Pb @, and 3.0 GeV on Pb [|29|], with 



standard CEM95 results (dashed histograms) and with calculations using modified pion absorption cross sections on 
quasideuterons, namely, W = 2 for the light C target, W = 3 for the medium mass Ni and Fe targets, and W = 4 
for the heavy Pb target (noted here as CEM98). 
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which are consistent with the macroscopic model used in the Moller-Nix mass calculations. We have 
significantly improved the treatment of rotational energy and have made a number of small refinements to 
the code like using actual updated values for the masses of elementary particles (for simplicity in CEM95, 
we used, e.g., m n = 140 MeV for pions of any charges and m n = m p = 940 MeV) and more precise 
values for physical constants. 

With these modifications and various coding improvements, we wrote a preliminary improved version 
of the CEM95 code called CEM97a, which is now incorporated in the MCNPX transport code ^Tt\\ . 

CEM97a is a preliminary release, since several improvements have been made after its creation and we 
continue the work in this direction, but already CEM97a provides a much better description of many data 
than the initial code, CEM95. As an example, in Fig. 12 we show excitation functions for production of all 
Xe isotopes for which we found experimental data from p+ 133 Cs interactions, part of our calculations for a 
medical isotope production study [|9|], using the standard CEM95 and with the improved version CEM97a. 
All available experimental data and calculations with the LAHET [|33|] code system (version 2.83) by 
K. A. Van Riper from [|] and with the HMS- ALICE code [TJ^] by M. B. Chadwick from the LA 150 
Activation Library [ |64] ] are shown as well, for comparison. One can see that CEM97a describes these data 
much better than CEM95 and LAHET, and the difference between predictions of the new and old versions 
of CEM increases as the isotopes produced become more neutron deficient, reaching a factor of 3 for 
120 Xe. Such a big difference between results from the new and old versions of the code is related mainly 
to using experimental nuclear masses and binding energies in CEM97a and the 40 year old approximations 
of Cameron [g] in CEM95. 

An example of the importance of using reliable values for nuclear masses can be seen also in Fig. 11. 
From this figure one can see that neutron spectra from Pb calculated with CEM95 agree very well with the 
data at low energies, in the evaporation parts of spectra. This is not necessarily good, since the calculations 
with CEM95 for Pb take into account competition between evaporation and fission, but the fission itself 
is not calculated; when a fission event occurs after the cascade and preequilibrium stages of a reaction, 
CEM95 stores this event to calculate a fission cross section, but then it stops further simulation of this 
event. That is, calculated spectra do not contain contributions of particles evaporated from fragments 
following fission. We are presently working to incorporate simulation of the fission processes themselves 
into the CEM, and we need some "room" for particles evaporated from fission fragments. As we see 
from Fig. 11, the new version (noted as CEM98) of the code provides us such "room"; calculated spectra 
agree very well with the data for intermediate and high energy parts of the spectra but lie below the 
measurements in the evaporation region by a factor of 2, leaving some place for neutrons evaporated from 
fission fragments. 

After creating CEM97a, we have further improved the code. We have improved the approximation to 
the level densities at the fission saddle point and have introduced other modifications, given separately in 
[]65t]. We also further modified the cascade modeling in the code by improving the treatment of photo- 
nuclear reactions and imposing energy and momentum conservation for the whole cascade stage of a 
reaction using experimental nuclear masses and binding energies for incident and emitted nucleons instead 
of the old approximations of CEM95. 

Our improvement efforts continue. We have written a subroutine for emission of complex particles at 
the cascade stage of reactions via the mechanism of their coalesence from emitted nucleons, have modified 
the preequilibrium and evaporation stages to allow emission of fragments with A > 4 from not too light 
excited nuclei, and are working now to develop an appropriate model for high energy fission. Among many 
other things, we need to improve the approximation of the inverse cross sections used in calculating particle 
emission widths and to treat more accurately a-emission at both the preequilibrium and evaporation stages. 
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Fig. 12. Excitation functions for the production of 127 120 Xe from p+ 133 Cs. Calculations with the new CEM97a 
code are shown by solid lines and with the standard CEM95 by short dashed lines. Results of calculations with the 



LAHET [33] code system (version 2.83) by K. A. Van Riper from Ref. \w are shown by long dashed lines and 



results of the HMS-ALICE code [39] calculated by M. B. Chadwick from the LA150 Activation Library [M] are 



shown by dotted lines. Experimental data are from the compilation [63]. 
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4. Summary 



In this paper, we have demonstrated the good overall predictive power of the modified version of the 
CEM as realized in the code CEM95. Then, to improve the agreement of its results with experimental data 
and to make it a better tool for applications, we modify it further, progressively incorporating features of 
previously neglected physics. So far, we have incorporated in the CEM new and better approximations 
for the elementary cross sections, imposed momentum-energy conservation for each simulated event, used 
more precise values for nuclear masses, Q-values, binding and pairing energies, corrected systematics 
for the level density parameters, studied and chosen the optimal approximations for the pion "binding 
energy", V n , for the cross sections of pion absorption on quasideuteron pairs inside a nucleus, for the 
effects of refractions and reflections, and for nuclear transparency of pions. We also make a number of 



refinements in calculation of the fission channel, described separately in Q65Q. 

As it was shown by a number of examples, the improvements to the CEM made so far clearly have 
increased its predictive power. Our work is not finished. Among improvements of the CEM which are 
of highest priority we consider development and incorporation of an appropriate model of high-energy 
fission, treating more accurately a-emission at both preequilibrium and evaporative stages, incorporation 
of a model of fragmentation of medium and heavy nuclei, possibly the Fermi breakup model for highly 
excited light nuclei, modeling evaporation and preequilibrium emission of fragments with A > 4, and 
improvement of the approximation for inverse cross sections. 

The problems discussed above are typical not only for the CEM, but also for all other similar models 
and codes, where they are also not yet solved. 
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